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L INTRODUCTION 


Large-scale motions of the atmospheres of planets and in the convection zones of rotating stars 
are strongly constrained by rotation, under the action of Coriolis forces, and gravity, which is 
manifest in the buoyancy forces that drive thermal circulations. The resulting flow structures 
are often surprising and continue to baffle scientists seeking fundamental understanding of such 
phenomena as the zonal bands of Jupiter, the origin of extremely high winds in the tropics and 
subtropics of Jupiter, Saturn and Neptune, the persistent differential rotation of the Sun, the 
complex patterns of convection in the slowly-rotating mantle of the Earth, and the rapidly 
rotating flows in the Earth's core that are thought to generate the Earth's magnetic field. 

Several idealized theories have been offered as explanations for some of these phenomena. 
Busse (1970, 1973, 1983) considers that Coriolis forces associated with planetary rotation 
dominate the large-scale (global) dynamics on the major planets. He hypothesized that motions 
in planetary atmospheres are constrained by the Taylor-Proudman theorem, which states that 
winds should be invariant along the axis of rotation for sufficiently high planetary rotation. His 
idea is sketched in figure 1. The axial invariance leads to a set of nested cylinders with 
convection occurring as columnar pipes inside the cylinders. Nonlinear interactions of 
individual convection columns can lead to large scale "zonal" winds (i.e. winds that are invariant 
with changes in longitude). The theories of Busse presume that the Taylor-Proudman theorem 
applies. However, this constraint is known to break down when stratification or turbulence (e.g. 
small-eddy viscosity) is strong enough. The precise nature of the breakdown of axial invariance 
due to turbulence is not known. A major goal of the USML2 experiments was to observe the 
parameter-space dependence of this process for comparison with scaling theory and 
computational simulations. 


2 THE PRINCIPLE OF THE GFFC 


The microgravity laboratory allows the study of fundamental physical processes that are 
masked or overpowered by gravitational influences in the terrestrial laboratory. One such 
physical force in dielectric liquids is that due to imposed electric fields. In a thermally stratified 
fluid, the dielectric polarization force is small compared to gravitational thermally induced 
buoyancy unless the field strengths and dielectric constant are unrealistically large (e.g. the 
electric field is well above the dielectric breakdown point). Therefore microgravity is required. 
The GFFC is an instrument based on electrohydrodynamic flow instability in a spherical annular 
geometry. It can be shown that for the particular electrodynamic system we use, the dominant 
flow instability is directly analogous to the buoyancy driven modes in rotating spherical shells. 
Thus observations of the evolution of the electrohydrodynamic instabilities, and their 
dependence on rotation rates and other parameters, show us how normal buoyant fluids should 
behave in the presence of heating and radial gravity. 

Detailed descriptions of the GFFC experiment have appeared previously (Hart et. al., 1986a, b, 
1990, 1991), so only a brief overview is given here. Fig. 2 shows a cross-section of the working 



cell. A dielectric silicone oil is contained in a hemispherical annular gap between a nickel 
sphere and a sapphire dome. The gap d = R o -R i is about 1 cm (see table 1). By using servo- 
controlled heaters imbedded in the inner sphere, and at the pole and equator of the outer sphere, 
various temperature distributions can be set up on the two surfaces bounding the working fluid. 
These are denoted T t (6) and T o (0) respectively, where 0 is the latitude. The simplest applied 
heating occurs when T a and T ( are constant with T i > T a . This latitudinally invariant 
convectively unstable heating distribution corresponds, for example, to astrophysically 
motivated simulations of stellar convection where nuclear burning at the core leads to a 
spherically symmetric destabilization of the convection layer. Latitudinally varying heating 
relates to planetary convection situations where the differential heating from the Sun either 
dominates or significantly affects the circulation in the fluid. Cases with a hotter pole or a hotter 
equator are both of interest. The latter has obvious connections to terrestrial planets like the 
Earth, but the geophysically relevant distribution of boundary temperature depends on how, in 
the case of planetary atmospheres, solar radiation is deposited with altitude. 

A large alternating voltage V a is imposed across the working fluid (a transparent but electrically 
conducting coating of Indium is applied to the inner surface of the sapphire dome for this 
purpose). The hemispherical annular gap shown in fig. 2 is continued below the equator to form 
a spherical capacitor. The large voltage (up to 10 kilovolts rms) gives rise to a substantial 
electric field in the fluid. There then occur both electromagnetic and gravitational buoyancy 
forces in the fluid which are given by 

F = ? £-l£’Ve + lvj^pfri-gP ,1) 

where q is the free-charge density, E(r,t) the instantaneous electric field which is a function of 
radius r and time t, p the fluid density, e the dielectric permittivity, and g the gravitational 
acceleration. E 2 is the dot product of the electric field vector. 

The use of an alternating potential with a frequency large compared to the charge relaxation 
time for the silicone oil working fluid leads to a negligible contribution from the first term on 
the right side of (1). Because the oil is incompressible and divergence free, the flows are 
governed by vorticity equations. Thus the dynamically relevant forces are those which survive 
after taking the curl of eqn. ( 1 ). Namely, 

VxF = p o aVTx V4> + £ o ^- VTx V£ 2 (2) 

where we have substituted the equations of state p = p o (\-aT) and e = e o (l -yT) for the 
density p and the permittivity e as a function of the temperature T. Here a is the coefficient of 
thermal expansion and y is the dependence of the dielectric constant on temperature. <p is the 
geopotential associated with whatever real gravity is present. Equation (2) shows that E 2 plays 
the same role as the geopotential in the dynamics. 
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Because E 2 = E*E has a time-average part that goes like r~* , V£ 2 is a radial vector decreasing 
like r -5 For the narrow gap used in the experiment the difference between r~ 2 and r 5 is 
relatively insignificant. However, because £ is quite small for non-polar fluids that can 
withstand a large applied potential, the electrodynamic buoyancy torques resulting from the curl 
of the polarization forces is small compared with terrestrial gravity. This smallness of the 
electrohydrodynamic polarization force is the main reason for placing the GFFC experiment in a 
microgravity environment. For example, table 1 gives typical values for the various 
parameters in the GFFC experiment. As shown in Hart et. al. 1986a, the effective 
electromagnetic gravity is 

„ =2e JW_ 

*“ ’O + Ppl'pr' 

where r is the radius and P is the ratio of the gap width d to the inner radius. Using a maximum 
voltage of 10 kilovolts and the values in table 1, it can be seen that the effective gravity is about 
0. 1 3 g's or less at the outer boundary, decreasing as V is lowered. 


3. PREVIOUS TERRESTRIAL EXPERIMENTS IN RELATION TO GFFC 

Previous laboratory studies of convection have been constrained to rotating layers between 
parallel plates oriented perpendicular to gravity (e.g. Rossby, 1969). This is because when 
rotating shells are used, terrestrial gravity is no longer everywhere perpendicular to spherical 
surfaces so that symmetry with respect to geopotential is lost and anomalous motions are 
generated. Even so, laboratory experiments in planar layers of rotating fluids have previously 
played a significant role in geophysical fluid dynamics. However, because terrestrial gravity and 
rotation are parallel in these experiments, only dynamics relevant to the polar regions of 
planetary atmospheres are addressed. Indeed it has long been a goal of fluid dynamicists to 
warp gravity in such a way as to allow a continuous variation of the angle between rotation and 
gravity vectors as is found in atmospheres and oceans of global extent. Prior to the GFFC this 
was only possible in a fluid of constant density, or in fluids with a very limited representation of 
thermal stratification. For example, Hart (1972) studied P-plane dynamics for two immiscible 
fluid layers in which the aforementioned angle varies from one side of the experiment to the 
other. 


4.1 CONTROL PARAMETERS AND EXPERIMENT PROCEDURES 

There are four fundamental dimensionless parameters for the GFFC. The Prandtl number (the 
ratio of kinematic viscosity to thermal diffusivity) is fixed at 8.4. The aspect ratio p is fixed at 
2.65. The Taylor number 


Ta = 4QV 4 / v 3 = 2.43x1 0 6 / r 2 

(where x is the rotation period of the cell) measure the effects of rotation relative to viscosity. 
The period ranged from 2 to 255 seconds. The Rayleigh number 
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Ra = gemY^Td* / kv =218.4V7y 2 , 


where VT is the applied radial temperature difference (degrees C) and V is the applied voltage 
in kilovolts, measures the relative effects of buoyancy and diffusion. MOST OF THE 
RESULTS IN THIS PAPER ARE GIVEN IN TERMS OF VOLTAGE AND PERIOD, SO 
THE ABOVE RELATIONS MUST BE USED TO CONVERT TO DIMENSIONLESS 
PARAMETERS WITH PHYSICAL SIGNIFICANCE. The temperature differences was 
typically about 15 degrees while the applied voltages ranged from 0 to 10 kV. In addition to 
these non-dimensional number relating the magnitudes of the experimental controls, there is 
also the prescribed distribution of temperatures on the inner and outer spheres. These scale with 
VT, but with spherically symmetric thermal boundary conditions no additional information is 
required to specify the system. On SL - 3 an extensive set of experiments with a latitudinal 
component to the heating were attempted. These were intended to mimic the fact that Jupiter’s 
atmosphere may be excited by both the Sun (giving a latitudinally varying heat input) and from 
internal heating (presumed to be isotropic). A few runs with latitudinal heating were attempted 
on USML-2 (we had hoped to do many). However, on USML-2 the thermal control system 
appeared not to be working properly, so such runs with large latitudinal gradients were curtailed. 
In any case a parameter 

H = pole to equator temperature difference / inner to outer radial temperature difference 
distinguishes between the various cases. In particular: 

H = 0 is called the “Solar” case, because it is thought that the solar atmosphere is not excited by 
a latitudinally varying heat source. 

H = 1 to 4 is called the “Jovian” case, for reasons stated above. Computational simulations and 
SL - 3 experiments were done for H = 1 to 3 . On USML-2 we had hoped to do runs 
with H = 3 and 4, but this proved impossible. 

Figure 3 shows a schematic of the GFFC instrument as flown on USML2. The cell is mounted 
on a turntable containing devices to regulate the temperature on the spherical shells. 
Visualization of the electrohydrodynamical buoyancy modes is obtained remotely, as any probes 
will perturb the desired spherical symmetry of the electric field and would find themselves in a 
very hostile electromagnetic environment. An aspheric lens assembly projects a Ronchi ruling 
onto the inner sphere. This is reflected back upon the ruled grid itself and the resulting image is 
then photographed on 16mm film. Upon traversing the working fluid, light rays are bent due to 
index of refraction fluctuations in the silicone oil that are in turn caused by temperature 
variations in latitude or longitude. Thus the film records images which depict radially averaged 
horizontal temperature gradients in the fluid. These reveal patterns of convection which arise 
under various external conditions: rotation rate Q, applied gravity (i.e. V), and heating 
distribution. During the filming, experiment performance data are recorded in digital (both 
decimal and binary) form on each frame. After the mission these data were read by a digital 
image acquisition system on the ground after the film is processed. The images of the 
convection cell were also recorded by a video camera. This video image includes four pieces of 
information concerning the experiment parameters. This is a very small fraction of the 
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parameter and performance data recorded on the 16mm film. Thus the film was the primary data 
storage medium for refined analysis after the flight, but the video system allowed for interactive 
experiments to be carried out during the USML-2 mission. 

Figure 3b illustrates the part of the GFFC hemisphere that is seen by the film and video camera. 
Although coverage includes the entire 90 span in latitude, only a limited range of longitude is 
covered. However, by combining 8 successive images that are taken each 45 degrees in 
longitude, we can present full hemispherical views. This works well for medium rotation rates 
(so the camera can take 8 images around a circumference quite rapidly), and for slowly evolving 
convection states (so that the convection is quasi-steady for one turn of the cell. The procedure 
for unwrapping the imagary is explained in APPENDIX 2. 

The generic experiment procedure was: 

a) Fix the external parameters of the experimental run. These are the rotation rate Q, the 
photography protocols, and the boundary temperature distributions. 

b) After an initial setup begin taking 16mm photographs through the backfocus Schlieren 
system. The frame rate was about 1 per rotation of the cell at high rotation rates (Q> 1 radian 
per sec.), making it possible to reconstruct a long time series of images to look for periodicities, 
aperiodicity, etc. 

c) Repeat with the new Ra (obtained by ramping or jumping the applied voltage V). 

d) After a range of Ra's are explored at fixed Q, Ti, and To, move on to cases with different 
rotation rates and heating distributions. Each Ra ramp experiment lasted 6 hours. In all, thirty 6 
hour runs were made. 


4.2 INSTRUMENT PERFORMANCE 


During the USML2 mission the downlinked astronaut monitoring information suggested there 
was a problem with the thermal controller. The actual temperatures obtained on the spheres did 
not match the commanded values very well. The binary LED data on the film, after using the 
postflight thermal calibration curves, confirmed this. Figures 4 and 5 illustrate the history of 
boundary temperatures for two 6 hour runs. Both runs have horizontally uniform commanded 
temperatures of 25 (outer sphere) and 40 (inner sphere) degrees. Because of calibration shifts 
(possibly during launch), the controller thought the outer sphere was too hot, and reacted by 
cooling it too much. However, during the run, as convection transferred heat from the inner to 
the outer sphere, the instrument cooler could not keep up and the inner sphere drifted upwards. 
Also notable is the non-uniformity of inner sphere temperature, approaching 10 percent of the 
applied gradient. APPENDIX 1 contains a summary of all the runs (which had film recording, 
see below). These detailed results are consistent with the comments made above. 

The astronaut monitors of the instrument panel suggested that experiments which commanded 
large imposed latitudinal temperature differences (with latitudinal variations as large or larger 


5 



than the requested radial variation) were in thermal error by order one amounts. In these cases, 
use of the postflight calibration indicated better performance than expected from the real-time 
monitoring data sent down during the flight, but still with considerable (30 percent) differences 
from the commands. This means that it will likely be necessary to use the detailed thermal data 
taken from the film to build useful interpretive computational models of the GFFC results 
because one cannot rely on the assumption that the commanded temperatures actually were 
attained during the USML2 runs. Unfortunately the film camera failed during run 19 (out of 
30), so that is impossible to get detailed quantitative results from the last 1 1 experiments 
because precise data on the thermal boundary conditions are not available for these runs. 
However, the downlinked video provides some qualitative information on flow states. 
Nonetheless, it was decided to concentrate our analysis and interpretation efforts on the 1 6mm 
film because of the much more complete instrument thermal records and much higher and more 
uniform temporal sequencing of the imagery. 


5 RESULTS FOR SPHERICALLY SYMMETRIC HEATING 

Runs 1 through 19 concentrated on the issue of the evolution of convection with spherically 
symmetric heating. Of these 19 runs, for which full 16mm film data were obtained, only runs 2 
and 19 had commanded non-symmetric (i.e. latitudinally varying) thermal boundary conditions. 

Figure 6a summarizes several of these spherically symmetric runs by showing individual 
snapshots of the convection on the parameter plane of the experiments. H = 0 denotes 
spherically symmetric heating. Each image shows the equivalent of an astronomer’s telescopic 
view of a star with the south pole at the bottom and the equator at the top (see figure 3b for a 
detailed discussion of the viewing). For high Taylor number (low rotation period, e.g. 2 
seconds) the convection initiates as prograde propagating polar disturbances, followed closely, 
as Ra is increased, by north-south oriented columnar (“banana cell”) modes (compare to figure 
1 ). 

The GFFC camera can be set to take a sequence of images each 45 degrees in longitude. If the 
motion is not evolving too fast in time, it is possible to take a set of 8 of these telescopic views, 
unwrap the peculiar birds-eye geometry, and merge the images into a full hemispheric (equal 
latitude line) projection. Figure 7 illustrates one example of this processing technique for a state 
with both propagating banana cell disturbances and non-periodic polar waves. A computational 
simulation of this state by Miller and Leslie is included for comparison. At weak to moderate 
supercriticality the comparison between simulation and experiment is reasonably good. 

As Ra moves higher so that the motions become turbulent, the polar convection eats northwards 
into the banana cell regime, basically eliminating the presence of axially symmetric elements in 
the motion. As the period of rotation is increased, the convection states become more 
complicated. The separation of polar and banana cells decreases and the eradication of banana 
cells with increasing Ra is more rapid. 

At very low Taylor number (long period of rotation) the convection is more transitory and 
patterns are less well defined. However, some runs showed evidence for persistent states with 
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relatively simple geometry. Figure 8 shows some examples from a run with a 48 second 
period. The voltage was commanded to ramp quickly to 1.5 kV, hold there, then run up to 3.5 
kV and back down down to hold at 1.5 kV again. During the first four hours of the run the 
states (fig. 8a, b) were irregular in time and space. However, the system finally settled down 
into a persistent wavenumber two propagating eddy (8c). One possible conclusion is that the 
system exhibits hysteresis (e.g. 8a and 8c are two fundamentally different states that are possible 
at essentially the same parameter setting). However, because the flows near 8a are not periodic, 
while 8c is, another interpretation is simply that it takes many hours for the fluid system to 
equilibrate and the images in 8a and 8b have not reached statistical equilibrium. 


5.1 COMPUTATIONAL SIMULATION OF SPHERICALLY SYMMETRIC FLOWS 

Figure 9 shows a regime diagram obtained from a 3D computational “Miller-Leslie” model (a 
variant of that described by Miller, et. al., 1992) using 30 radial points, 62 latitudinal points and 
45 longitudinal wavenumbers. The model boundary conditions correspond to the experiment, 
except that purely symmetric temperatures (e.g. averaging fig. 4 in latitude) were used. At high 
Taylor number the states are similar with prograde propagating polar modes coming in first as 
Ra in increased, and clear indication of banana cells in the mid-latitudes and tropics. The states 
at low Taylor number appear to be more stable (less transitory) than in the experiments. 
Multiple states (or hysteretic behavior) were found in the simulations. For example there are 
two qualitatively different but persistent states of motion at Ta= 1000, Ra = 20,000. In general, 
the simulations using this model do a good job at describing the motions at high Ta and low to 
moderate Ra, but seem to give overly simple states at low Ta. At low Ta the onset of 
convection occurs at a lower Ra, and the supercriticality, at a fixed voltage (3 kV, for example), 
is much higher so that there may be resolution problems in some of the low Ta runs. Also, at 
this time it is not known how sensitive the computed flows are to modest latitudinal variations in 
the boundary temperatures. 


5.2 THE ONSET OF CONVECTIVE INSTABILITY 

A new camera action protocol, implemented on the USML2 flight that gives a closely spaced 
time series of images at a fixed longitude permits an accurate determination of the onset of 
convection (i.e. the fundamental horizontal symmetry breaking instability). Figure 10 shows the 
observed onset points. As illustrated in figures 6 and 7, there are two types of modes, prograde 
polar and prograde propagating banana cells, the latter having strong north-south alignment. 
The experimental curves represent fits through data from 9 runs at different rotation rates. At 
high Ta polar modes come in first as Ra is raised, while at lower Ta the banana cells onset first. 
Figure 7b shows the onset wavenumber and propagation rates for banana cells. The 
wavenumber drops rapidly as the rotation period increases, as does the prograde propagation 
rates of the banana cells. 

Onset data can be compared with predictions of computational simulation, and this provides a 
useful validation of the latter because the weak motions associated with the instability tests most 
of the terms in the model while not requiring enormous 3D resolution. Figure 10 shows that the 
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Miller-Leslie model comes quite close the experimental result for the 4 second period, although 
the separation between polar and banana modes is smaller than observed and the onset points are 
slightly higher. More simulations that include a fit to the actual (not smoothed) thermal 
boundary conditions are planned. The Glatzmaier (1983, see also Hart, Glatzmaier, Toomre, 
1986) model has higher resolution, but because it is a pseudo-spectral code, it does not 
implement the rigid boundary condition of the GFFC present at the equator (see figure 2). The 
tips of the banana cells should generate Ekman layers long the equatorial barrier, and should be 
more damped in the presence of this wall. Thus it is not too surprising that the Glatzmaier 
model predicts that banana cells should arise first, at a lower Ra. 


5.3 TRANSITION LATITUDES FOR TURBULENT FLOWS 


A very important result from the GFFC experiments, obtained by combining data from the 1 8 
runs done with symmetric heating, is an estimate of the parametric conditions required to 
observe axial (banana cell) convection in highly nonlinear and turbulent flow. Such results 
cannot be obtained computationally over the entire range of the GFFC parameter values. What 
is the latitidunal limit for banana (axially aligned) structures in the flows? If this critical latitude 
shrinks to zero, then the object will not have a banana cell regime. If the results from the 
experiments conform to an asymptotic scaling relation, then it might be possible to use the 
GFFC results to predict whether or not a rotating star or planet will have aligned giant cells in 
its atmosphere. Fig. 10 shows data for one rotation rate. As the voltage increases the 

banana cell domain shrinks towards the equator. For V > 8 kilovolts, the basic planetary 
rotation is too weak to align the motion in the north-south direction. 

Attempts to derive a scaling in this complex geometry typically involve some perhaps tenuous 
assumptions, but the GFFC results can be used as an independent check. One simple idea is that 
the transition boundary is related to the existence of more constrained geostrophic convection 
equatorward of a transition latitude 0 , and we take 0 to be the latitude at which the convection 
has unit Rossby number. The Rossby number measures the relative strengths of advection and 
Coriolis forces and is given by R a =U/ 20. sin(0)d where U is a characteristic convection 
veloc ity. If we take the convection velocity to be estimated by buoyancy-induced free fall, then 
= yjg em yATd and setting the Rossby number equal to 1 gives the scaling relation (for fixed 
Prandtl number): 

sin(0) R~ V2 T a m « (r V)-'. (3) 


An alternative model, using estimates of U from the numerical P-convection model of Brummell 
and Hart (1992) that includes local boundary curvature in the dynamics, gives a somewhat 
different scaling: 


sin (0)~R a - VA T a 


15/24 


y-V2 r -5/4 


( 4 ) 


In the second result the transition latitude falls faster with V because in the Brummel-Hart free 
slip P-convection model, the motions develop high speed jets by mean flow instability. The 
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zonal flows dominate the convection and the high velocities associated with these jets lead to a 
prediction that unit Rossby number is reached sooner . 

Fig. 1 1 shows that the data are consistent with the simple thermal Rossby number scaling. The 
data are not consistent with the P-convection scaling. This is in concert with the observation 
from GFFC that no high speed jets (which would be associated with anomalous disturbance 
propagation rates) were observed. Also, the P-convection model predicts vacillatory states 
involved with mean flow instability. Such periodic pattern pulsation states were not seen in the 
USML2 GFFC results either. Thus full sphericity seems to detract from the potential for strong 
jet formation. This result needs to be tempered by the possibility that the GFFC experiments, 
conducted with a Prandtl number of 8.4, may not accurately predict what might happen in a low 
Prandtl flow, or in atmospheres with forcing parameters well outside the GFFC experimental 
range. If, however, the Prandtl number dependence turns out to be weak (as checked, perhaps, 
by computational simulation), the scaling (3) may allow determination of the flow regimes to be 
expected in rotating stars and planets where the applied temperature difference (or 
supercriticality) is known. Unfortunately whereas (3) is independent of diffusion parameters (a 
plus!), if one must use observations of emitted heat flux in order to estimate the Rayleigh 
number (through a Nusselt number scaling), this then re-introduces the thermal diffusivity into 
the scaling. Future work in the application of such scaling relations as (3) and (4) to 
astrophysical situations will determine whether or not the obstacle of having only a poor 
knowledge of small-eddy diffusivity is significant. 


6 CONCLUSIONS 


The GFFC experiments on USML2 returned an extensive data set on the nature of 
convection in rotating spherical shells with nearly horizontally symmetric applied thermal 
boundary conditions. The main results extracted to date have many facets. Some results are 
discouraging with respect to our initial expectations, but some results were surprising and should 
provide fuel for more theoretical investigation and computational simulation that will increase 
our general understanding of rotating convection in spherical shells: 

On the down side: 

1) There was lack of experimental evidence for vacillatory states of convection that had been 
predicted in some simple theoretical models. This may be a geometry and Prandtl number 
problem. 

2) There was no robust evidence for banding or axisymmetric jet formation as hoped for in 
order to explain differential rotation on the Sun or the multiple zonal jets on the giant planets. 
We did observe differential propagation of patterns (prograde , retrograde, prograde, moving 
from equator to pole), suggesting three different zones of instability. We do not know whether 
or not mean zonal motion plays a role in the differential propagation (perhaps computational 
simulations could help here), but whatever mean flows are present seem weak compared with 
the overturning motions of the convection eddies (i.e. the propagation rates are slow compared 
with the pattern fluctuation timescales). 
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These experiments seem to rule out rapidly rotating banana cell convection as a cause for these 
phenomena when operating over the parameters range studied here. Thus, if banana cell 
convection is the root cause, it must be operating at a significantly different setting. A 
possibility is that parameters may be such as to place the system in the banana republic regime 
(see figure 1 1), but at a much higher Rayleigh and Taylor number. Experiments could address 
this situation only in a geometry with a much larger gap width. 

3) Although SL-3 returned some very interesting data on H unequal to zero flows (that were 
successfully interpreted using theory and numerical models), USML-2 experiments had 
difficulty in replicating or extending the SL-3 experiments because of thermal control 
problems. 

On the plus side: 

1) There is good agreement between computational simulations of the onset of instability and 
the experiments. 

2) The observation of much more complicated and transitory convection at low rotation than 
found in the computational simulation models suggests a sensitive dependence on physical 
boundary conditions, or on computational resolution. An unraveling of this situation by the 
modelers will no doubt improve our ability to simulate flows of this type. 

3) Unanticipated patterns of convection at modest rotation rates (i.e. 16 second and 48 second 
periods) also will prove to be a useful challenge to computational codes. The mixture of quasi- 
axisymmetric motions at mid latitudes with non-axi symmetric motions at the pole and the 
equator is particularly surprising. The degree to which this may be related to the presence of 
the equatorial baffle in the GFFC instrument remains to be sorted out. 

4) Approximate verification of a simple scaling law for the existence of columnar convective 
turbulence in highly stressed shells. This relation and its simple physical interpretation should 
be useful in categorizing the convection state of rotating geophysical objects like planets and 
stars. 

All in all, the two flights of the GFFC have provided much data on rotating convection over 
about 80% of the parameter range in the original design. There is still much work necessary to 
interpret physically the results from USML-2, particularly in the low rotation cases. 


7 SUGGESTIONS FOR FURTHER ANALYSTS 

Future analyses and comparisons with models should focus on how weak latitudinal variation in 
boundary temperatures may influence the dynamics, whether or not simulations with higher 
resolution can recover some of the low rotation GFFC results, and on further characterizations 
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of the important high Rayleigh number turbulent states observed in rotating spherical shells. Some 
of the more interesting questions to pursue include: 

1) Can the transition latitudes for banana convection be reproduced computationally? 

2) What caused the bizarre patterns of convection at low rotation (e.g. appendix 2, figs 4 -7)? 

3) Were 6 hour runs long enough to reach equilibrium? The models currently suggest much faster 
adjustment times than were observed experimentally. Are the models internally incorrect in their 
numerics, etc., or are they missing some essential physics? Perhaps the convection patterns at low 
rotation are very sensitive to the precise distribution of temperature on the boundaries, not only 
spatially, but temporally. In the USML-2 GFFC runs the temperatures did drift with time (by of 
order 10-20 %). Was this enough to keep some of the experiments from settling into a steady 
pattern? 

On the experimental side, the use of electrostatically generated radial gravity in the configuration 
of the GFFC has proved its benefit and feasibility. However, to extend the range of parameters into 
more practically significant ranges, with Taylor and Rayleigh number thousands of times larger 
than that possible in the GFFC, an entirely new concept will be needed. This parameter range will 
likely elude modelers for many years, but is of much interest geophysically. We still don’t under- 
stand the manner in which large scale coherent banana cell structures may re-emerge in a flow 
whose transition latitude has long ago (in a parameter sense) reached the equator. The intriguing 
possibility that sufficiently turbulent convection with many small eddies may organize itself into 
large scale structures (by collective action of the small scale eddies) might be addressed by labora- 
tory experiments with a much deeper (x 20 or more) gap. The use of ferromagnetic fluids merits 
study in this regard. 
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Table 1. Experiment parameters and characteristics of Dow Coming 0.65 centistoke silicone oil. 


Nominal Fluid Properties (Dow Corning 0.65 cs 200 Fluid) 


Ambient density 

P 

760 kg/m 3 

Expansivity 

CL 

1.34 xlO -3 °C _I 

Kinematic viscosity 

V 

6.5 x 10 -7 m 2 /s 

Thermal diffusivity 

K 

7.7 x 10“ 8 m 2 /s 

Ambient permittivity 

e 

2.5e 0 

Vacuum permittivity 

^0 

8.90 x 10 -12 farad/ 

Dielectric variability 

7 

1.29 x 10" 3 °C _1 

Conductivity 

a 

10 -12 mho/m 

Heat capacity 

c 

1.7 x 10 3 joules/kg 

Dissipation loss factor 

<P 

«4x 10 -5 

Nominal Experiment Parameters 

Rotation rate 

Q 

0-3 rad/s 

Radial temp, difference 

at; 

0-25 °C 

Voltage (r.m.s.) 

V 

0-10 kV 

Voltage frequency 

Cl) 

300 Hz 

Inner radius 

Ri 

2.402 cm 

Outer radius 

K 

3.300 cm 

Gap 

d 

0.908 cm 

Aspect ratio 

P = RJd 

2.65 

Prandtl number 

Pr 

8.4 


Table 1 
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TABLE 2: GFFC (6 HOUR) RUNS PERFORMED ON USML-2, 1995 

Run Scenario Rotation CorruTianded / Obtained voltage Film Camera 


Number 

Number 

Period 

TC6 

TIE 

TOP 

TIP 


Maqazine 

Action 

Comments 

1 

1 5 

2 

25 

40 

25 

40 

2->l0 

1248 

3 

Solar Ramp 




27.4— >33.6 

40 

27.9->34.3 

41 





2 

1 3 8 

48— ►S 

20 

30 

30 

40 

3 

1240 

2 





29.8 

39.7 

30.3 

40.7 







8 

20 

30->25 

30->35 

40 

5— >10 


3 





33.6— >33.8 

43-5 — ►3'8.7 

34.3->36.0 

44.2— >40.7 





3 

2 2 

4 

25 

30 

40 

45 

2— >10 

1244 

3 

Baro Ramp 




32.2-* 34 .2 

3 1.1 — »32 

37.5 ->40.6 

43.6->45.2 




H-3 not very well 

4 

1 7 

4 

25 

40 

25 

40 

2->10 

1244 

3 

Solar Ramp 




25.6— >31 

39.9— >40.1 

263->31.1 

40.7— >40.1 





5. 

1 1 2 

4 

25 

40 

25 

40 

1 — >9 

1244 

8 

Solar Ramp 




27.2->30.8 

40 

27.8->32.4 

40 





6 

1 8 2 

8 

25 

40 

25 

40 

0->9 

1245 

3 

Solar Ramp 




25.6— >31.0 

40 

262->31.0 

40.7 





7 

1 8 1 

1 6 

25 

40 

25 

40 

0->9 

1245 

2 

Solar Ramp 




27.9— >31.0 

40 

27.9->32.6 

405 





8 

1 7 1 

3 

25 

40 

25 

40 

l->9 

1246 

9 

Solar Ramp 




25.6— >31.0 

40 

26— >32.4 

42.4-f40.7 





9 

1 7 2 

6 

25 

40 

25 

40 

l->9 

1246 

9 

Solar Ramp 




27.8— >31.1 

40 

27.9->32.6 

40.7 





1 0 

1 9 8 

48 

25 

40 

25 

40 

052->8 

1246 

2 

Mantle Ramp 




25.6— >30.0 

40 

26.3 — >31-1 

40.7 





1 1 

1 9 1 

96 

25 

40 

25 

40 

0->9 

1335 

2 

Mantle Ramp 




25.6— >29.8 

40 

262->31.1 

40.7 





1 2 

1 4 5 

4 

25 

40 

25 

40 

0->9 

1335 

9 

Solar Ramp 




25.6— >31.1 

40 

263->375 

42.4— >40.7 





1 3 

1 9 2 

128 

25 

40 

25 

40 

052— >5 

1335 

1 

Mantle Ramp 




27-5— >29.4 

40 

27.9->30.4 

40.7 





1 4 

1 9 3 

240 

25 

40 

25 

40 

052->5 

1335 

1 

Mantle Ramp 




25.8— >29.8 

402— >40.0 

263->30.8 

40.7 





15 

1 94 

40 

25 

40 

25 

40 

0.8->35 

1335 

2 

took for Mantle Hyst. 




26.2— >29.5 

39.8-H40.0 

26.8->29.9 

40.7 

->0.8 



01.6kV 

1 6 

1 9 5 

128 

25 

40 

25 

40 

0.8->35 

062 

1 

Look for Mantle Hyst. 




29— >29.8 

40->39.8 

29.6— >30.4 

40.7— >42.4 

->0.8 



01.6kV 




->29.2 


->29.6 

—►40.7 





1 7 

1 9 6 

48 

25 

40 

25 

40 

3,1,2, 

062 

2 

Look for Mantle Hyst. 




29.8— >30.4 

402-»39.8 

30.6— *3 1.1 

40.7 

1.3,2 




1 8 

1 9 7 

128 

25 

40 

25 

40 

2,3, 1, 

962 

1 

Look for Mantle Hysl. 




29.4— >30.8 

39.*_>40.i 

302->31.1 

424->40.7 

2, 3,2 





—►42.4 
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\jn 

Scenario 

Rotation 

Commanded / Obtained 


Voltage 

Film 

Camara 


mber 

Number 

Period 

TOE 

TIE 

TOP 

TIP 


Magazine 

Action 

Comments 

1 9 

1 6 1 

48 

30 

45 

45 

55 

0.44— *4 

962 

2 

Jovian Lo V Ramp 




36.4-07.5 

45.8-444.8 

43.4—4455 

535-06 




■No film last 1/3? 







-*55.3 





20 

2 9 

48 

25 

35 

35 

45 

2— *10 


2 

Jovian Hi V Ramp 




27.9-01.1 

35.2 

36— *34.3 

45.6— *43.9 




No film 







-445.6 





2 1 

1 6 3 

2 

30 

45 

45 

55 

3— *9 

1247 

8 

Jovian Hi V Ramp 




35.9-07.5 

452— *44.8 

42.8-445.6 

552-035 




No film last 1/6 







-052 





22 

1 6 4 

8 

30 

40 

50 

55 

2— *9 

1244A 

8 

Baro Ramp, H-3 




35.8-+402 

39-»4Q.7 

495— *503 

535 




No film last 273? 





—>4051 

-*47.1 




r 


23 

1 6 5 

4 

30 

45 

45 

55 

3— *10 

1245A 

Jr 

Jovian Ramp 




35.6—07.5 

42-445.0 

455-*43.9 

552-4535 




Film for 5 min. 

24 

1 5 7 

2 

25 

40 

25 

40 

2 — ► 10 


8 

Solar Ramp 




26-01.1 

40 

26.8—02.4 

42.4-440.7 




No film 

25 

1 5 8 

2 

25 

40 

25 

40 

2-410 

1335A 

9 

Solar Ramp 




26.6-00.8 

40 

275—431.1 

42.4 




Film for 20 min? 

26 

1 7 3 

6 

25 

40 

25 

40 

1— *9 


8 

Solar Ramp 




25.6-01.1 

40 

263-02.7 

42.4-440.7 




No film, Video? 

27 

1 6 9 

4 

30 

40 

50 

55 

2-*9 

1246A 

8 

Baro Ramp, H-3 




37-440.2 

402 

48.6— *503 

535-4552 




Film for 5 min? 






—*49.1 






28 

2 2 0 

1 6 

30 

40 

50 

55 

1-49 


2 

Baro Ramp, H-3 




38.3->395 

403 

47.9-*503 

—*47.1—4495 

552-4532 




No film 

29 

2 1 9 

128 

25 

40 

25 

40 

2, 3.2 


1 

Look for Mantle Hyst. 




25.6-427.9 

40 

26.3-*27.9 

42.4 




Q 2 kV, No film 
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TABLE 3. GFFC RUNS PERFORMED ON SPACELAB - 3, 1985 




(commanded parameters shown beiow) 




Scenario 

Rotation 

TOE 

TIE 

TOP TIP 


High 

Comments 

Number 

Period 

(degreeC) 

(degreeC) 

(degreeC) (degreeC) Voltage Ramp 



(sec) 





(kV) 


35 

2 

25 

35 

25 

35 

3.6 -> 10 

6 hour Solar Ramp H = 0 

40 

48 

25 

40 

25 

40 

5.0 -> 10 

3 hour Solar Ramp H = 0 

16 

2 

25 

30 

40 

45 

0-> 10 

3 hour "Jupiter 1 * w/ latitudinal gradient H = 3 

20 

2 

25 

30 

45 

50 

5 -> 10 

3 hour "Jupiter" w/ latitudinal gradient H = 4 

21 

2 

25 

30 

45 

50 

1 ->3 

3 hour "Jupiter w / latitudinal gradient H = 4 

23 

16 

25 

30 

40 

45 

1 -> 3 

3 hour "Jupiter w / latitudinal gradient H = 3 

9 

2 

25 

35 

25 

35 

0-> 1.6 

5 hour solar with dye injection, H = 0 

5 

250 

25 

35 

25 

35 

0-> 1.0 

5 hour basic instability run, H = 0 

12 

2 

30 

25 

45 

40 

0->10 

3 hour ocean, stably stratified 

13 

2 

30 

25 

45 

40 

0 -> 2 

3 hour ocean, stably stratified 

60 

20 

38 

35 

38 

35 

5.44 fixed 

20 minute spin up run, stably stratified H = 0 

43 

2 

30 

45 

30 

45 

5 - 6 - 5 - 6 

voltage up down up jumps, H = 0 

44 

2 

30 

40 

30 

40 

2 - 3 - 2 - 3 

voltage up down up jumps, H = 0 

24 

16 

25 

35 

35 

45 

5 -> 10 

3 hour "Jupiter w / latitudinal gradient H = 1 

25 

16 

25 

35 

35 

45 

1 -> 3 

3 hour "Jupiter w/ latitudinal gradient H = 1 

41 

48 

25 

35 

35 

45 

5 -> 10 

3 hour "Jupiter w f latitudinal gradient H = 1 

42 

48 

25 

35 

35 

45 

1 -> 3 

3 hour "Jupiter w / latitudinal gradient H = 1 

18 

2 

25 

35 

35 

45 

5 -> 10 

3 hour "Jupiter w/ latitudinal gradient H = 1 

22 

16 

25 

30 

40 

45 

5 -> 10 

3 hour "Jupiter w / latitudinal gradient H = 3 

16 

2 

25 

30 

40 

45 

0-> 10 

3 hour "Jupiter" w / latitudinal gradient H = 3 

17 

2 

25 

30 

40 

45 

0 -> 2 

3 hour "Jupiter w / latitudinal gradient H = 3 

19 

2 

25 

35 

35 

45 

1 -> 3 

3 hour "Jupiter" w / latitudinal gradient H = 1 

20 

2 

25 

30 

45 

50 

5 -> 10 

3 hour "Jupiter w / latitudinal gradient H = 4 

21 

2 

25 

30 

45 

50 

1 -> 3 

3 hour "Jupiter w/ latitudinal gradient H = 4 

34 

2 

25 

35 

25 

35 

1.4 -> 3.2 

6 hour solar ramp, H = 0 

38 

16 

25 

35 

25 

35 

3.6 -> 10 

6 hour solar ramp, H = 0 

4 

20 

25 -> 40 

25 

25 -> 40 

25 

0 

Basic calibration, no applied voltage 

5 

250 

25 

35 

25 

35 

0-> 1.0 

5 hour basic instability run, H = 0 

46 

2 

25 

45 

25 

45 

5 -> 10 

3 hour solar, H = 0 

37 

16 

25 

35 

25 

35 

1.4 -> 3.2 

6 hour solar, H = 0 

39 

48 

25 

40 

25 

40 

0.48 -> 2.0 

3 hour solar, H = 0 
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Figure 2. Schematic cross-section of the Geophysical Fluid Flow Cell. The gap contains 0.065 
centistoke silicone fluid. 


OJ 





Cross-section of the optical imaging system. 


Figure 3a. Instrument cross-section. The optics provides a Schlieren visualization of convection 
in the spherical convection cell. A beam splitter diverts the image to both a 16mm and a video 
camera. 
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Pole 



Fig. 3b Sketch of the view of the GFFC sphere that is seen by the optical system. Much 
as a telescope views a planet the viewing area extends from the equator to the pole 


19 




GFFC Scenario 145 Outer Sphere Temps. 



Figure 4. Thermal history at the walls for scenario 145 (4 second basic rotation period, high voltage 
ramping from 0 to 10 kiloVolts). The upper panel shows the inner sphere values at the latitudes shown, 
while the lower panel illustrates the behavior of the outer sphere at the latitudes shown. 
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GFFC Scenario 194 Inner Sphere Temps. 



Exp. Elapsed Time in Hours 


GFFC Scenario 194 Outer Sphere Temps. 



Exp. Elapsed Time in Hours 


Figure 5 Thermal history at the walls for scenario 194 (48 second basic rotation period, high voltage 
ramping from 0 to 3 kV and back down again). The upper panel shows the inner sphere values at the 
latitudes shown, while the lower panel illustrates the behavior of the outer sphere at the latitudes 
shown. 
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Rotation Period (sec) USML2 
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Ta = 2.43x1 0 6 /period 2 

Figure 6. Regime diagram showing qualitative results for several Rayleigh number ramp up experiments with symmetric heating. 
Each image is a single snapshot with the south pole at the bottom and the equator at the top. The solid curve shows the onset of 
banana cell, the dashed the onset of Dolar waw convection F^rh imoo^ ^n-penonHe tn rw uo — i — i— 





Figure 7a. Unwrapped convection planforms for scenario 145 (4 second rotation period, symmetric 
heating, 2.16 kilovolts). The top two panels show motions as viewed by E-W fringes (left) that are 
sensitive to longitudinal temperature gradients, and N-S fringes (right) that are sensitive to latitudinal 
gradients. The bottom panel shows radially averaged temperatures form the Miller-Leslie model. 
Note the general agreement on structure and equivalence of disturbance longitudinal wavenumber. 
In all frames the view is from pole (center) to equator (limb) with equal spacing of latitude lines. 
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Figure 7b. Propagation characteristics near onset for spherically symmetric 
heating and the basic rotation rates (periods) as shown. Top panel shows 
the near— critical azimuthal wavenumber . Bottom panel shows near-critical 
(prograde) wave propagation speed in radians per second. 






Figure 8. Unwrapped views (as in figure 7) of scenario 1 94 (48 second rotation, symmetric heating). 
N-S fringes 9 sensitive longitudinal thermal gradients) for values of voltage shown (1.46, 3.16, 
1.5 kV). The first two images (top left and right) are during the first three hour ramp up, while the 
bottom image (the only stable pattern) is at about hour five. The first image is during the first 
voltage hold period, the second during the ramp up to 3.5 kV, and the third is just after the end of the 
second voltage hold period. 
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Ta = 2.43 x 10 6 /period 2 

Figure 9. Regime diagram showing computed planforms in the Miller-Leslie model (Fourier decompostion in longitude, grid 
points in the meridional plane). Boundary conditions, r-dependent gravity, centrifugal acceleration, and geometry as in GFFC. 


Onset Rayleigh Number vs. Taylor Number 


1 00 Rotation Period (sec) 1 0 



Figure 10 . Critical curves from the GFFC experiments. Polar modes become the most unstable 
for Taylor numbers in excess of about 20,000. The results from computational simulations are 
shown as points on the figure 
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APPENDIX 1 

DATA REDUCTION SYSTEM 

TEMPERATURE AND HIGH VOLTAGE RECORDS FROM 
USML-2 MISSION 16mm FILMS 
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GFFC Data Reduction 


The process of extracting data from the GFFC experiment film, and producing information 
easily read and viewed using a PC, required three networked computers, a stale of the art high resolution 
video camera, 2 video capture boards, a laser videodisc recorder, a computer controlled film projector, 
along with device control and image analysis software developed at the University of Colorado. 

USML2 and Spacelab 3 film of the Geophysical Fluid Flow Cell was converted to video on a 
modified stop action movie projector. Computer control of the movie projector allowed for accurate frame 
positioning and illumination level control of the experiment film to be viewed with a Kodak Megaplus 
model 4.2 black and white video camera. The Kodak video camera is stale of the art equipment with a 
2K X 2K square pixel aspect ratio and an anlog ouput signal with ten bits ( 1024 gray levels) of dynamic 
range. Film illumination control was required because of intensity variations between picture type 1 and 
picture type 2 images. Capture of the very high resolution video signal from the Kodak video camera was 
achieved with a Coreco Oculus 500 video frame grabber for PC compatible systems. The video capture 
card served two purposes. First the card was used to reformat the video signal to an RS-170 three 
component (RGB) video signal compatible with a Sony laser videodisc recorder (LVR-3000N). Each 
image of the 135,000 plus film frames was stored on laser videodisc. Second the film annotation for the 
USML2 mission was digitized and stored on computer hard disk for processing later. All elements of 
image acquisition, projector control .data compression , and data storage were controlled by a C program 
written just for the GFFC image recording system. 

After digitizing, compressing ,and storing USML2 annotation to computer hard disk the image 
files were passed via an Ethernet connection to an SGI Infinite Reality Power Onyx for conversion from 
image to ASCII data. Even with compression ratios of 10: 1, the image data alone occupied over 10 
gigabytes of disk space. To convert the experiment image annotation data to a more readable ASCII 
format, a program was developed using the IDL graphics language. This language was selected for it ease 
of use with image data. Annotation image quality variations required a program that would adapt to light 
level changes, image shifts, distortion, background noise, and misbehaving experiment LEDs. The 
interactive nature of the IDL programming environment made it easy to test program modules for proper 
operation. This reduced program development time and increase our confidence in accurate character 
recognition. 

Annotation data came in two formats. Seven segment LED numeric characters were recognized 
by determining the on/off state of the segments. This proved to be faster, more adaptive, and easier than 
matching the numeric characters to templates (0 to 9) Dot Matrix data also was deciphered in a similar 
manner by checking the on/off state of the LED's and converting the binary data to an ASCII floating 
point value. Difficulties in the image recognition process occurred because of image variations across a 
frame, from frame to frame, and from scenario to scenario Unfortunately the program was not always 
able to recognize the annotation. For this reason, the output data file contains an error code which 
indicates the confidence of the data recognition. 

The complete and calibrated data set for the USML2 mission was converted from 10 gigabytes of 
compressed image data to approximate 10 megabytes of ASCII data. In addition to converting the image 
files to ASCII format, the IDL program also corrected the data using post flight calibration curv es The 
curves were developed by placing the experiment apparatus in a temperature controlled chamber and 
recording measured thermistor sensor outputs vs. actual values for the range of temperatures the 
experiment was operated at during the USML2 mission. Similar calibration curves were produced for the 
hemisphere voltage and turntable rotation period by comparing measured or command values to actual 
values. 


In order to make viewing of the GFFC experiment images easy’ and accessible to anyone with a 
PC, the images were converted from laser videodisc to computer readable Video for Windows format. 
Using a Miro video DC30 capture card the GFFC experiment images on laser videodisc were digitized 
and written in the Motion Jpeg file format. Because this format is hardware specific, the Motion Jpeg 
scenario movie files were converted again to a more standard Video for Windows format using Adobe 
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Premiere 4.0. In addition, Adobe Premiere was used to add titles and remove the unwanted picture type 3 
images from Spacelab 3 data. The 320 X 240 X 16bit .AVI movie files produced were compressed using 
the 'CRAM 1 compressor, resulting in a data set from both USML2 and Spacelab 3, which could be fit on 
six compact disks . Also included with the image data is the ASCII experiment annotation data set. 

Using the windows base aviplay.exe program both the scenario images and calibrated experiment 
parameter data can be viewed simultaneously. 
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Thermistor Positions on Hemispherical Shells 


October 2, 1996 




A = Heat Flux, Watts 
B = Tj e ,°C 
C = T| p , °C 


D = T oe ,°C 
E = T 0 p, °C 
F = HV, KVolts 


G = Time, Min. Sec 

H = Scenario Frame N 
K = Turntable Rotatio 
Period, Sec 
L = Picture Type 
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GFFC Data Formats 


• Files in /data directories 

Data is in ASCII format and is readable using most word processors (e.g. Notepad). 

IDL format descriptor used to write the data 
format = 14 TR3 15(F5.2,TR2) 13 TR4 F5.3 

Format description: 

Four character integer, three spaces, fifteen floating point values with two spaces separating each value 
( Each number is five characters wide including the decimal point, and two characters follow the decimal 
point), three character integer, four spaces, five character floating point value with three characters 
following the decimal point. 

All files in /data directories contain 120 lines of data. 

Column labels for data file 

FRAME #, TOP, TOE, TIP, TIE, HV, IS 15, 1S25, IS35, IS45, IS55, IS65, IS75, OS25, OS45, OS65, RPRD, ETIME 


• Files in /usml2 directory on disk 1 

Data is in ASCII format and is readable using a word processor capable of loading large documents 
(e.g. Word). 

IDL format descriptor used to write the data 
format = 14 TR3 15(F5.2,TR2) 13 TR4 II 


Format description: 

Four character integer, tliree spaces, fifteen floating point values with two spaces separating each v alue 
( Each number is five characters wide including the decimal point, and tw'o characters follow' the decimal 
point), tliree character integer, four spaces, 1 character integer. 


Column labels for data file 

FRAME #, TOP, TOE, TIP, TIE, HV, IS 15, IS25, IS35, IS45, IS55, IS65, IS75, OS25, OS45, OS65, RPRD, ECODE 


FRAME# 

TOP 

TOE 

TIP 

TIE 

IS15 

IS25 

IS35 

IS45 

IS55 

IS65 

IS75 

OS25 

OS45 

OS65 

RPRD 

ETIME 

ECODE 


= USML2 film frame number (frame numbers wrap around at 9999) 

= Temperature at the outer pole 
= Temperature at the outer equator 
= Temperature at the inner pole 
= Temperature at the inner equator 
= Inner Hemisphere Temperature at 15 degrees 
= Inner Hemisphere Temperature at 25 degrees 
= Inner Hemisphere Temperature at 35 degrees 
= Inner Hemisphere Temperature at 45 degrees 
= Inner Hemisphere Temperature at 55 degrees 
= Inner Hemisphere Temperature at 65 degrees 
= Inner Hemisphere Temperature at 75 degrees 
= Outer Hemisphere Temperature at 25 degrees 
= Outer Hemisphere Temperature at 45 degrees 
= Outer Hemisphere Temperature at 65 degrees 
= Commanded Turntable Rotation Period in seconds 
= Experiment Elapsed Time in hours 

= Error Code ( 0=data OK, 1= film annotation overlap-data suspect. 3= data corrupt) 
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• Files in /movies directories 

Video for Windows compatible 
Frame size = 320 X 240 
Depth is 16bits 

Compressor "CRAM" , Microsoft Video 1 [32] 
Keyframe every frame 
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THE FOLLOWING FIGURES SHOW THERMAL AND VOLTAGE HISTORIES FOR ALL USML-2 

EXPERIMENTS WITH COMPLETE l6mm FILM RECORDS 
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GFFC Scenario 164 Hemisphere Temps. 
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GFFC Scenario 181 Hemisphere Temps. 
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GFFC Scenario 182 Hemisphere Temps. 
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GFFC Scenario 191 Hemisphere Temps. 
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APPENDIX 2 


UNWRAPPING USML-2 IMAGES ONTO A FULL 

HEMISPHERE 
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APPENDIX - 2: GENERATION OF 360 DEGREE LONGITUDE VISUALIZATIONS 

FROM THE GFFC IMAGERY. 

An outline of the mapping routine used to generate the full hemispherical flowfields: 

For certain camera actions it is possible to extract a sequence of eight images separated 
longitudinally by 45° on the spherical convection cell. If the temporal evolution of the 
observed flow structures is not too rapid during the sampling period for this sequence, 
then the images may be used to produce a composite image depicting the flowfield across 
the entire hemisphere as seen from the pole. 

The developed image mapping scheme consists of three stages. In the first stage each 
recorded image is projected back onto the surface of the sphere, thereby defining the flow 
structures as a function of longitude and latitude at that radial location. This process 
utilises the nature of the wrap around optics used to record the images. These optics 
cause the paths of all the observed radially propagating rays to become refracted on 
emerging from the convection cell, so that they are parallel to the optical axis of the 
camera. Then in the second stage, the above results are used to map each image in order 
to attain its view as seen from the pole. Since orthographic projections cause an 
excessive amount of image distortion near the equator, a mapping routine that generated 
polar plots in which radial distance was proportional to latitude was implemented An 
example of one such transformation is given in Fig. 1 b. Lastly, in the third stage the 
eight mapped images within the sequence are overlaid to produce the final composite. In 
order to produce a smooth transition between images at successive longitudinal locations a 
weighted mask was applied to each The effect of this mask is to smoothly vary the 
intensity of the image azimuthally from a maximum at its centre to zero at a position of 
± 45° . The weighting function was chosen such that the total weight of all the individual 
contributions at any location was equal to 1.0, thus ensuring a uniform illumination across 
the composite. An example of one such masked image together with the final resulting 
composite are shown in Fig. 1 c+d respectively 

The operation of the above scheme was validated by observing its effect on the lines of 
equal latitude etched onto the surface of the glass convection cell. As can be seen from 
the results, especially in cases highlighting north to south intensity gradients, these marks 
are faithfully mapped onto circles of corresponding latitudes 

Figures 2 - 7, following, illustrate full-view maps obtained for several of the runs in 
table 2 of the summary section. 
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Figure 1 a) Recorded image; b) Mapped image; c. Masked image; d, Final Composite. 
60 





,***S*£JP 




W 


mr 








5/i023 
V: 2.8( 


& -4* 


fel-V j 




JT*- 


1 5/i02647— 1 
kV: 3.04000 


1 5/i02648-2 
kV: 3.04000 


I ^£2*5^ 


■&■*, E 




r 

I 


Sv 


1 5/i02807-l 
kV: 3.13000 


1 5/i02808-2 
kV: 3.13000 


Figure 2. Unwrapped views of convection in scenario 1 5 with values of the voltage as indicated. 
Left columns are N-S fringes (sensitive to E-W gradienst) and the right columns show E-W 
fringe setting (sensitive to N-S gradients). Note the development of banana cell aligned with the 
rotation axis, following onset of polar disturbances. The applied kilovoltage is shown for each 
pair of images. 
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Time:26 min Time:26 min 

kV:2. 16000 kV:2. 16000 


Figure 3. Unwrapped views of convection in scenario 17 with values of the voltage as indicated. 
Left columns are N-S fringes (sensitive to E-W gradients) and the right columns show E-W 
fringe setting (sensitive to N-S gradients). Convection is similar to figure 2 , except note that the 
banana wavenumber is lower. 
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Figure 4. Unwrapped views of convection in scenario 181 with values of the voltage as 
indicated. Left columns are N-S fringes (sensitive to E-W gradients) and the right columns 
show E-W fringe setting (sensitive to N-S gradients). Although the rotation period is long 
compared with figs. 3 and 2 (16 seconds, vs. 4 and 2 respectively), banana cells still form, 
although the polar disturbances are more pronounced. As the voltage is increased the polar 
modes remain beautifully organized and stronger than the banana cells. 
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kV: 1.00000 kV: 1.00000 



182/i 14546-1 1,11 \W^' 182/i 14547-2 

kV: 1.04000 kV: 1.04000 



Figure 5. Unwrapped views of convection in scenario 182 (8 second period) with values of the 
voltage as indicated. Left columns are N-S fringes (sensitive to E-W gradients) and the right 
columns show E-W fringe setting (sensitive to N-S gradients). The response is qualitatively 
similar to figure 4. 
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182/i 159 13-1 
kV: 6.12000 




182/115914—2 
kV: 6.12000 




kV: 3.13000 kV: 3.16000 



Figure 6. Unwrapped views of convection in scenario 194 (48 second period) with values of the 
voltage as indicated. Left columns are N-S fringes (sensitive to E-W gradients) and the right 
columns show E-W fringe setting (sensitive to N-S gradients). Note the mixture of 
axisymmetric convection and NS striations that evolves into strong steady wavenumber 2 mode 
as kV is raised. 
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APPENDIX 3 


DETAILED COMPARISON OF GFFC RESULTS WITH 
VARIOUS COMPUTATIONAL SIMULATIONS 
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Comparisons to numerical simulations: 


Two numerical models were used to make some general predictions of GFFC behavior 
before the USML2 flight and then to compare with actual GFFC results after the flight. The 
models are set up with the GFFC spherical geometry, electro-gravity, and the actual fluid 
parameters. Both models attempt to be direct numerical simulations of the experiment in that they 
have no additional dissipation beyond Laplacian molecular diffusion. A DNS model must have 
sufficient spatial resolution (i.e. a large enough number of grid points or spectral modes) so that 
nothing except molecular diffusion processes occur between the grid points in finite-difference 
formulation or at larger than the highest modeled wavenumber in spectral formulation. In 
addition, smaller spatial resolution implies smaller temporal resolution for numerical stability. 
Hence the more turbulent the flow, with small high-speed flow structures, the more expensive the 
computation. As shown below, the numerical model results are most like the GFFC results near 
convective onset as expected. 

The model by Miller and Leslie (Miller, et. al., 1992) is 3-dimensional finite-difference in 
latitude and radius, and spectral in longitude. The post-flight runs discussed below were run with 
30 radial points on a stretched grid that gives more points in the boundary layers, 60 latitudinal 
points and 45 longitudinal waves. The boundary conditions are able to mimic the temperature 
boundary conditions of the GFFC as well as the no-slip insulating boundary at the equator. The 
model by Glatzmaier (Glatzmaier, 1983 and Hart, et al. 1986) is spectral in latitude and longitude 
with a slightly better spatial resolution, and Chebyshev-spectral in radius. The spectral technique 
is numerically more accurate but does not implement the no-slip equatorial boundary of the GFFC 
apparatus. Instead, Glatzmaier requires symmetric solutions between northern and southern 
hemisphere, equivalent to a free-slip velocity boundary condition there. 

Post-flight comparisons to individual scenarios: 

Example 1) Solar case, 4 second rotation period, no intended latitudinal temperature gradient 

The voltage and temperature inputs to the numerical model were obtained from the GFFC 
experiment temperature and voltage readouts recorded on the 16 mm film (scenarios 145, 1 12, 
and 17). These conditions were simplified slightly (for numerical convenience) from the actual 
conditions of the USML-2 runs (see Appendix 1 for details on the actual temperature 
distributions) to: 

Tinner = 38. 5C 
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Touter = 18C from 0 to 70 minutes, then a linear amp 18 - 24C over 70 - 200 minutes, and 
finally a linear ramp 24 - 26C during 200 - 360 minutes. 

Linear voltage ramp from 0.806 - 4.95 kV during 0 - 230 minutes, then a ramp of 4.95 - 
9.50 kV during 230 - 360 minutes. 


Experiment: 

In this scenario as time progressed the forcing voltage increased in the GFFC experiment, 

and convection began at voltage 1.57 kV near the pole with a longitudinal wavenumber of 
about 17. This pattern traveled prograde at about .0035 radians per second. When the 
voltage increased at 1.85 kV we get the typical north-south aligned “banana cells”. These 

banana cells of wavenumber —17 traveled prograde at speed ?. At slightly larger voltage the 
banana cells lengthened toward the pole and interacted with the polar convection. Further 
increase in voltage saw the polar convection region extend equatorward and erode the banana cell 
pattern. At the same time these equatorial cells became nonlinear in amplitude, decreased in 
average wavenumber, had energy in a wider range of wavenumbers and traveled more slowly. 
Examples of the observed patterns for increasing high voltage values are shown in Figure 1 The 
unwrapping technique was used to reconstruct the spatial dependence of the full hemispheric 
pattern from successive images taken every 1/8 revolution. In la the polar convection has begun 
and the equatorial bananas are just starting in the upper-left quadrant. By lc the banana cells are 
interacting with the polar convection. 

Miller-Leslie Model: 

Two example images from the simulation are shown in Figure 2. In both, the temperature 
is shown on a color scale which runs from blue to red to green to yellow from cold to hot. Recall 
that the GFFC images show either the north/south or east/west temperature gradient rather than 
temperature directly as in the model. The patterns are barotropic with depth, D. Figure 2a is near 
onset of the model banana cells and 2b has nonlinear bananas interacting with the polar 
convection. In Figure 3 time-longitude plots of temperature at mid-radius at latitudes of 66 and 
1 5 degrees show the onset values and phase speeds of the polar and equatorial convection 
respectively. The polar convection begins in the model at 2. 10 kV, significantly higher than 
observed in the experiment, with a prograde-traveling wavenumber 5 pattern. Within a few 
minutes this simple pattern breaks into a complicated set of convection rolls which move both the 
simple pattern becomes . At 2.12 kV the banana cell convection begins, closer to the observed 
onset. The pattern is a mixture of wavenumbers 18 and 19, also close to the experimental values. 
The prograde longitudinal phase speed starts rapid (~9.3 deg./min. at time = 75 min. and voltage 

= 2 16 kV) and as the amplitude grows and wavenumber shrinks, slows to ~5.3 deg./min. at time 
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= 85 min. and voltage = 2.34 kV). In summary, the model calculates the onset of polar 
convection at too large a voltage, does better with the equatorial convection onset, finds the onset 
banana cell wavenumber very well, and does ? with the phase speed of the banana cells. The 
decrease in wavenumber and phase speed of the banana cells for increasing high voltage is 
consistent with the GFFC observations. 

Glatzmaier model: 

This model calculates the onset at of equatorial convection at about 1.6 kV, before the 
polar convection begins. The free-slip boundary at the equator in this model presumably damps 
the convection there less. This onset voltage for bananas is close to but lower than that observed 
in the experiment but completely misses the initial polar convection mode. Figure 4 shows 
contour plots of temperature and vertical velocity at D=0.5 for this model. The onset 
wavenumber for the bananas, 18, is consistent with both GFFC and the other model. As the 
voltage increases to the range of 3 to 5 kV, this fully spectral model does a qualitatively better job 
than the M-L model, predicting more GFFC-like, higher horizontal wavenumber structures. 

Example 2) Mantle convection case, looking for possible hysteresis, 48 sec rotation period, no 
intended latitudinal temperature gradient, high voltage ramping up, holding, up again, down, 
holding again at same value, and then down again. 

The voltage and temperature inputs to the numerical model were obtained from the GFFC 
experiment temperature and voltage readouts recorded on the 16 mm film for scenario 194. 

These conditions were simplified slightly for numerical convenience to: 

Tinner = 39C 

Touter = 18C from 0 to 70 minutes, then a linear amp 18 - 23C over 70 - 180 minutes, and 
finally a linear ramp 23 - 21C during 180 - 360 minutes. 

Linear voltage ramp from 0.66 - 1.47 kV during 0-70 minutes, 
constant voltage of 1 .47 kV during 70- 106 minutes, 

linear voltage ramp from 1.47 to 3.49 kV during 106 - 175 minutes, 
linear voltage ramp from 3.49 to 1.47 kV during 175 - 246 minutes, 
constant voltage of 1 .47 kV during 246 - 28 1 minutes, 

linear voltage ramp from 1.47 to 0.64 kV during 281 - 360 minutes. 

Experiment: 

Preflight numerical simulations with the M-L model had predicted hysteresis could be 
observed in the slowly-rotating cases. Indeed, in the first half of this scenario when the voltage 
was increasing then held fixed and then increasing, structures like Figures 5a, d were observed 
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with long rolls which extended a large fraction of the pole-to-equator distance and a circumpolar 
ring at high latitudes. For larger voltage (above about 2.5 kV) higher wavenumber, shorter rolls 
were observed near the pole along with a low-latitude (~20 deg.) ring which had smaller scale 
rolls attached to it (Figures 5b, e). For decreasing voltage a very stable wavenumber 2 ring was 
observed (Figures 5c, f) rather than the high-latitude ring and north/south rolls. This wavenumber 
2 structure propagated prograde. It is not clear from the GFFC data that the pattern during the 
two constant voltage times was truly statistically steady and thus hysteretic or whether the entire 
run was a very long transient. Certainly the flow pattern was qualitatively different for increasing 
or decreasing high voltage. 

Miller-Leslie Model: 

The convection began right away in the simulation just as in the scenario 194 experiment, 
with an unsteady meandering of convection rolls that weren't obviously affected by the slow 
rotation. During the first constant voltage section, the pattern settled into a high-latitude ring 
with three long north/south rolls attached to it (Figure 6a). This pattern was stable and 
propagated. It is similar to the GFFC pattern of Figure 5a, d but more regular and more stable 
At 2.0 kV the number of spokes declined to two and then at 2. 1 kV to one. From 2.4 kV up to 
the maximum and then back down to the 2nd holding voltage at 1 .46 kV the convection pattern 
looked like Figure 6b with a cool ring at mid-latitude and hot ring at high latitude and small-scale 
structures on each. Both rings and the waves on them are time-dependent. This pattern is 
reminiscent of Figures 5b, e but differs in the position of the rings and in the smallness of the other 
rolls. Furthermore, after less than 10 modeled minutes in the 2nd holding voltage, the entire 
convection pattern dissipated into axisymmetric flow rather than the stable 2 lobe structure of 
Figures 5c, f. In the GFFC experiment, the flow was convecting until the end of the run Thus 
the modeled pattern is different for the two holding voltages, but not at all like the experiment In 
general, the separate patterns are also more stable with less noise in the model and the transitions 
between patterns are shorter than in the experiment, so that separate flow regimes or qualitative 
patterns are more easily identified in the model. 

Glatzmaier Model: 

Qualitatively, the flow regimes are similar to GFFC during the first 2/3 of the run while the 
voltage increased, held steady, increased, and then decreased (Figures 7a-d) Again this model 
predicted axisymmetric flow before it was observed in the experiment, but not until the last hour 
of the scenario. The robust 2-lobed structure that dominated the 2nd half of the GFFC run was 
not calculated by this model Figures 7e), f). 
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Example 3) Mantle convection case, looking for possible hysteresis, 128 sec rotation period, 
no intended latitudinal temperature gradient, high voltage ramping up, holding, up again, down, 
holding again at same value, and then down again. 

The voltage and temperature inputs to the numerical model were obtained from the GFFC 
experiment temperature and voltage readouts recorded on the 16 mm film for scenario 195. 

These conditions were simplified slightly for numerical convenience to: 

Tinner = 39C 

Touter = linear ramp from 18 - 23C over 0-180 minutes, and then a linear ramp 23 - 22C 
during 180 - 360 minutes. 

Linear voltage ramp from 0.66 - 1.46 kV during 0-53 minutes, 
constant voltage of 1.46 kV during 53 - 105 minutes, 

linear voltage ramp from 1.46 to 3.46 kV during 105 - 172 minutes, 

linear voltage ramp from 3.46 to 1.46 kV during 172 - 230 minutes, 

constant voltage of 1 .46 kV during 230 - 279 minutes, 

linear voltage ramp from 1 .46 to 0.64 kV during 279 - 360 minutes. 

Experiment: 

Again, pre-flight numerical simulations had predicted that hysteresis might be observed in 
this slowly-rotating case. The patterns during the first half of the run until the end of the long 
voltage hold section (Figures 8a, g) are similar to the patterns observed in Scenario 194, long, 
unsteady, convection rolls which span most of the sphere from pole to equator. During the 
continuing voltage ramp from 1.46 to 3.5 kV (Figures 8b, c and 8h, i) the pattern switched to a 
high latitude ring with 2 spokes going all the way to the equator, very much like the pattern in the 
scenario 194 simulation by M-L (Figure 6a). Further into the scenario, the pattern became a 
soccer-ball of long convection rolls which persisted until the end of the run (Figures 8d-f, j-1). 
Within this basic pattern the larger voltage cases had more and sharper "noise" rolls. Like 
scenario 194, the convection pattern differed at like voltages for increasing and decreasing high 
voltage. Again it is not clear that either pattern was statistically steady or just stuck in a long 
transient. 

Miller-Leslie Model: 

The convection began right away in the simulation just as in the scenario 195 experiment, 
with an unsteady meandering of polar convection rolls that weren't obviously affected by the slow 
rotation (Figure 9a). During the first constant voltage section, the modeled flow became sub- 
critical and the pattern died away to axisymmetric flow (Figure 9b). When the convection began 
again at ~ 2.7 kV, the pattern of Figure 9c was maintained until early in the 2nd voltage holding 
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level at 1 .46 kV when the convection died away to axisymmetric flow again. This simulation is 
quite different from the GFFC runs. Qualitatively, the modeled scenario 194 looks most like 
GFFC scenario 195. 

Glatzmaier Model: 

Again the flow regimes are similar to GFFC very early in the run. This model also 
predicts very weak convection and nearly axisymmetric flow beginning during the 2nd voltage 
holding section (Figure lOe). 


Pre- and Post-flight Regime diagrams: 

Before the USML2 flight, images from the Space Lab 3 GFFC experiment and both 
numerical models were combined in a j$lot in (Rayleigh number, Taylor number) space in order to 
have a basic idea of what flows we should expect. Figure 1 1 shows a post-flight version of such a 
diagram. The USML2 GFFC data is quite extensive for H = 0 (see table 2, lead section) and only 
3 of the runs are shown. Comparable diagrams are shown in Figures 12-20 for SL3 data, the M-L 
model and the Glatzmaier model. 
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kV: 1.93000 



kV: 2.08000 kV: 2.08000 



kV: 2.16000 



kV: 2.16000 


Figure 1. Unwrapping from 4 second rotation period, solar case, no latitudinal temperature 
gradient. Scenario 17 (same as 145, but with camera action 3). Grating in N/S orientation (a-c) 
and E/W orientation (d-f) . a,d) Polar convection with equatorial banana cells just beginning to 
grow, b, e) Banana cells with longitudinal wavenumber ~17. The pattern travels prograde (CW). 
Banana cells with longitudinal wavenumber ~16. 
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Figure 2) Miller-Leslie simulation of Scenario 145. D is the nondimensional height in the fluid layer of the 
temperature field shown. "Average" is the vertically averaged temperature field. This polar projection has latitude 
mapped linearly to radius as in the unwrapped GFFC images a) 214 kV b) 2.41 kV. 
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Figure 3) Miller-Leslie simulation of Scenario 145. Temperature contours for all longitudes as a function of 
time at mid-height and latitude a) 66 deg. and b) 15 deg. 

Temperature at D=0.5 and 66 deg. Latitude 
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Figure 6) Miller-Leslie simulation of Scenario 194. D is the nondimensional height in the fluid layer of the 
temperature field shown. "Average" is the vertically averaged temperature field. This polar projection has latitude 
mapped linearly to radius as in the unwrapped GFFC images, a) 1.64 kV during the increasing voltage segment b) 
1 .65 kV during the decreasing voltage segment. 
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figure 7) Glatzmaier simulation of Scenario 194. Polar 
of contours of temperature at mid-height at the end of each hour 
1.6 kV, e) 1.2 kV, f) 0 6 kV. 


projection (with latitude linearly mapped to radius) 
during the run. a) 1 4 kV, b) 1 .9 kV, c) 3 2 kV, d) 
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Fieure 9) Miller-Leslie simulation of Scenario 195. D is the nondimensional height in the fluid layer of the 
temperature field shown. "Average" is the vertically averaged temperature field. This polar projection has latitude 
mapped linearly to radius as in the unwrapped GFFC images, a) 0.75 kV veiy early in the ran b) L93 ^kV ust 
after the holding voltage, c) 2.97 kV. This basic pattern remained until the convection died away at 1 .46 kV during 


the 2nd holding voltage. 
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(c) 
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Ta = 2.43x1 0 6 /period 2 

Figure 11. Regime diagram from the USML2 GFFC flight experiment showing qualitative results for several Rayleigh number ramp 
up experiments with symmetric heating (H= 0). Each image is a single snapshot with the south pole at the bottom and the equator at the 
top. The solid (dashed) curve shows the onset ol banana cell (polar wavy) convection. Each image corresponds to the dot to its left, in 
vertical order. 
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Figure 12-19. Pre USML2 flight regime diagrams organized by heating function, H (or organized by “SL3”, “M-L model”, 
“Glatzmaier model?”) 
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Figure 13. Computational simulations: Miller model H=0. 
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Figure 14. Computational simulations: Glatzmaier model H=0. 
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Figure 15. Summary of SL-3, H=1 states. 
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Figure 16. Miller Model H = 


Glatzmaier Model H= 

Rotation Period (sec) No Centrifugal Force 
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Figure 17. Glatzmaier Model H = 


Rotation Period (sec) Miller Model H=3 
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Figure 18. Miller Model H = 




Rotation Period (sec) Glatzmaier Model H=4 
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Figure 19. Glatzmaier Model, H = 
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